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Based on the slave-boson theory of the two-dimensional 
t ~ t' — J model, we calculate the superconducting conden- 
sation energy for optimally doped and overdoped high Tc 
cuprates at finite temperatures using a renormalized random 
phase approximation. The contributions come from the mean- 
field part and the antiferromagnetic spin fluctuations. In the 
presence of neutron resonance peak at (7r,7r), the latter is 
shown to have similar temperature and doping dependences 
as the difference in antiferromagnetic (AF) exchange energy 
between the normal and the superconducting state. This dif- 
ference has been proposed to be related to the superconduct- 
ing condensation energy by Scalapino and White. The total 
condensation energy, however, is about 1/2 smaller than the 
proposed AF exchange energy difference and shows a more 
rapid decrease as the temperature rises. The doping depen- 
dence of the condensation energy is found to be consistent 
with experiments. In particular, near zero temperature, our 
result shows a good quantitative agreement with experiments. 

PACS number: 74.25.-q,74.25.Bt,75.40.Gb,74.72.-h 



I. INTRODUCTION 

One of the important aspects to understand the mech- 
anism of superconductivity is to investigate the origin of 
the condensation energy, which is the free energy differ- 
ence between the normal and the superconducting(SC) 
states. In conventional BCS superconductors, the con- 
densation energy is accounted for by the change in the 
ion kinetic energy between the normal and the SC phases, 
as shown by Chester Q. In high-Tc superconductors, it 
is now widely accepted that the strong electronic interac- 
tions rather than the electron-phonon interactions may 
be responsible for the superconductivity. However, no 
consensus has been achieved on the detail pairing mech- 
anism, consequently on the origin of the condensation 
energy [|-|. 



In analogy with the phonon-mediated electron pair- 
ing, Scalapino and White Q proposed that if the pairing 
is mediated by the antiferromagnetic (AF) exchange in- 
teraction, the condensation energy will be proportional 
to the change in the AF exchange energy between the 
normal and the SC states. This gives a direct connec- 
tion between the condensation energy and the dynam- 
ical spin susceptibility x (or the spin structure factor 
5(q,c^) = Imx(q,^)/[l + exp(-w/T)]) as ff, 

Q poo 1 

En-Es^PJ2 / — [inix'^(q,^)-inix^(q,^)] 



x[cos((7^) -|-cos(qj,)]- — -. 

1 — exp(—u/I ) 



(1) 



The dynamical spin susceptibility can be probed directly 
by neutron scattering experiments. The most prominent 
feature in these experiments when the temperature de- 
creases below the SC transition is the appearance of a 
sharp resonance peak at the AF wavevector Q — (tt, tt) 
(tt resonance) in the SC state |@-^. From this obser- 
vation, Demler and Zhang further argued that the 
appearance of the neutron resonance will cause the an- 
tiferromagnetic exchange energy to be lower than that 
in the normal state. Without resorting to any micro- 
scopic models, they estimated the exchange energy to be 
around 0.016J based on the experimental data for the 
spin structure factor and the g-width of the resonance 
peak. The estimated condensation energy is roughly of 
the same magnitude as that obtained by experiments |Tc| ] 
and therefore accounts for a large part of the conden- 
sation energy. The above estimate is based on Eq.(|^) 
without including the kinetic energy. Furthermore, it 
is limited to zero temperature and with the assumption 
that the spin susceptibility in the normal state has neg- 
ligible contribution to Eq.(l). Its validity thus relies on 
further tests against experimental data. In particular, as 
pointed out at the end of Ref. p| , a useful test of this idea 



1 



is to compare the theoretical estimate with experiments 
at finite temperatures. Experimentally, because the net 
difference of the AF exchange energy is small, it will re- 
quire extremely careful measurements to check Eq.(p. 
Only recently, a study of the temperature dependence is 
performed. In Ref. |ll] , the neutron scattering measure- 
ment shows a similar temperature dependence between 
the condensation energy from the AF exchange interac- 
tions and the resonance peak intensity for optimal doped 
YBa2Cu306-i-£c. In view of the experimental situation, it 
would be useful to carry out calculations at finite temper- 
atures based on an existing high- Tc model. Recently, the 
t — t' — J model has been shown to produce the resonance 
peak hence it would be interesting to verify its 

prediction for the condensation energy. 

In this paper, we examine the temperature and dop- 
ing dependences of the condensation energy in a mi- 
croscopic model - the two-dimensional (2D) t — t' — J 
model. To avoid the complications due to the appearance 
of the pseudo-gap in the normal state, we shall confine 
our calculations in the optimal and overdoped regions 
iQ. Theoretically, the condensation energy can be cal- 
culated from the difference in the thermodynamic poten- 
tial between the normal and the SC states because the 
difference in free energy in the SC state is equal to the 
difference in thermodynamical potential JTst . In order 
to include the tt resonance peak, we calculate the ther- 
modynamical potential in the linked cluster expansion 
via a renormalized random-phase approximation (RPA) 
as defined in Ref. This approach has been previ- 

ously tested in different contexts [ p^p^Jl^ , p7[ | and has 
been shown to capture many important features of the tt 
resonance peak. Here we shall further investigate how it 
would predict for the temperature dependence of the con- 
densation energy. We show that the condensation energy 
coming from the AF exchange correction is nearly the 
same as Eq.(0) due to the presence of the neutron reso- 
nance peak. Furthermore, the condensation energy from 
the AF exchange correction has the similar temerature 
dependence as the resonance peak intensity for both op- 
timally doped and overdoped systems. This is consistent 
with the experimental results observed in the optimally 
doped sample However, the total condensation en- 
ergy after including the mean-field free energy difference 
between the normal and SC states exhibits a more rapid 
drop in its temperature dependence than that calculated 
with Eq.(|]). We also investigate the doping dependence 
of the total condensation energy and find it in agree- 
ment with the experiment. In view of the appearance 
of the pseudo-gap in the normal state for underdoped 
cuprates ||l8| where a part of the condensation energy has 
gained above Tc this result indicates that the con- 
densation energy for the optimally doped and overdoped 
cuprates should show more rapid temperature drop than 
that for the underdoped cuprates. 

This paper is organized as follows. In Sec. II, we in- 



troduce the model and define notations. The condensa- 
tion energy is derived using the linked cluster expansion. 
In Sec. HI, we calculate the temperature and doping de- 
pendences of the condensation energy and compare them 
with that calculated based on Eq.(|^). Finally, we give a 
concluding remark in Sec. IV. 



II. THERMODYNAMIC POTENTIAL AND 
CONDENSATION ENERGY 

We start with the 2D t — t' — J model which reads. 



j^CjCT h.c. 



<ij>,(7 <ij>' ,cr 

-fJ^S,.S„ (2) 



<ij> 



where, < ij > denotes the nearest-neighbor (n.n.) bond, 
< ij >' the next n.n. bond and Si = ^cl^fa/gCi/g. In 
the slave-boson method, the physical electron operators 
Cia- are expressed by slave bosons bt carrying the charge 
and fermions ficr representing the spin; Cia = bffic- In 
the SC state, we consider the order parameters 
fi]fjl ~ /u/jT with the c?-wave symmetry and 

= J2a < fiaf]<y >= Xo, in which bosons condense 
hi — >< hi >— y/S (S is the hole concentration) Then, 
the mean-field Hamiltonian in the SC state is, 

i?™ - E ^'^/La- - E + ^-^o 

ka k 

+2NJ'{xl + Al), (3) 



where = —2{St + J' Xo)[cos{kx) + cos(fcy)] — 

ASt' cos{kx) cos{ky) — fi is the dispersion for fermions, and 
Afe = 2J'Ao[cos(fc^) - cosiky)], with J' = 3J/8. In the 
optimal and overdoped regions, the mean-field Hamilto- 
nian for the normal state is obtained by setting the SC 
gap Aq = |2^. The mean- field parameters xo, Aq 
(Ao=0 in the normal state) and the chemical potential fi 
for different doping S are obtained from a self-consistent 
calculation p9| ]. 

It has been shown |T2|jr^ , |20| | that the spin suscepti- 
bility at the mean-field level represented by the first 
fermionic bubble in Fig.l fails to describe some important 
features observed by inelastic neutron scattering experi- 
ments, such as the resonance peak The underlying 
reason is because the mean-field theory underestimates 
the spin fluctuation near (tt, tt) ||21[]. Therefore, it is nec- 
essary to go beyond the mean field approximation. For- 
mally, this can be done by perturbing around the mean- 
field Hamiltonian, i.e., we write H = Hm + H' , and treat 
H' as a perturbation. In principle, all the fluctuations 
are included. However, different selection of subset dia- 
grams may result in different kinds of fluctuations. For 
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the spin fluctuation, the usual random phase approxima- 
tion selects a series of ring diagrams as shown in Fig.l. 
The resulting spin susceptibility can be resummed as 



(q, = 



(4) 



where 7^ = cos{qx) + cos{qy), xo is the unperturbed spin 
susceptibility which comes from the fcrmionic bubble. 
The parameter rj is formally introduced to keep track 
of the renormalization of vertex and its value is one in 
the usual RPA approach. However, the spin suscepti- 
bility calculated using Eq.(^ with i] = 1 leads to an 
AF instability for doping 5 < 6c ^ 0.22, which is much 
larger than the experimental data around Sc — 0.02. It 
indicates that there exists other fluctuations which acts 
as to suppress this overestimation. In the renormalized 
RPA approach [|2|, one models the suppression of 77 by 
treating 77 as a phenomenological parameter, whose value 
is determined by setting the AF instability at the ob- 
served value S = 0.02. For the material parameters we 
use, we found = 0.34 for the SC state and rj = 0.15 
for the normal state. Using 77 = 0.34 for the SC state, 
Brinckmann and Lee have explained the resonance peak 
and the incommensuration structure observed in neutron 
scatterings, and some of the present authors have quan- 
titatively accounted for the doping dependences of the 
peak/dip/hump structure in the angle-resolved photoe- 
mission spectra [ p^ and the tunneling spectra in a NS 
junction 

The above considerations can be carried over to the 
calculation of the thermodynamic potential. We consider 
the same set of ring diagrams in the linked cluster expan- 
sion of the thermodynamic potential as shown in Fig. 2, 
and find it is given by p2[ 



l + I^T]JjgXo {l,l^n) 



Here, we have used the renormalized strength 77 J for the 
vertex. The parameter v is used to keep track of the 
number of times that the potential appears in the per- 
turbation expansion. Since v enters in the same way as 
a coupling constant, the integration can be taken outside 
of the summation over icun- Performing the summation 
over iujn, we get 



An 



with 



^E 



dv / — Imx^ (7/,q,w)7,coth( — , 
Jo ^ 



X+ {v, q, uj) = 



xt 



(q^t^) 



l + vriJlqXt (q^^) 



(6) 



(7) 



In obtaining Eq.(0), we have made use of the symmetry 
Imx(q,-w) = -Imx(q,w). 

The condensation energy H'^v/^tt per unit cell is equal 
to the difference of the thermodynamic potential(free en- 
ergy) between the normal and the SC states. In our 
approach, it includes the contribution from the AF ex- 
change correction AJl' and that from the mean-field 
Hamiltonian. The free energy at the mean-field level in 
the SC state can be easily calculated from Eq.(||) as. 



--yin(l 
+2NJ\Al + xl) 



(8) 



with Ek = Ve| 



Its counterpart in the normal 



state Fq can be obtained from Eq.(8) by setting Ag — 0. 
To include the contribution from x, we note that the 
mean-field free energy already contains the Hartree-Fock 
contribution, i.e., the one-loop term in Fig. 2 j23j. There- 
fore, we should substract it to avoid double counting and 
arrive at, 

HJ.V 



Stt 



^ V / dv I — [Imx^(i^,q,w) - Imx'^(i',q,w)] 
^ g Jo Jo 



X7, CO 



X7,coth(— ) + Fo~-Fo^ 



(9) 



where, we have abbreviated x^ ^-s X- The first term 
of the right hand side of Eq.(9) comes from the AF fluc- 
tuation correction as shown in Fig. 2, while the second 
term is the Hartree-Fock contribution which should be 
substracted from the first term. Finally, — Fq is the 
free energy difference at the mean-field level. 

Before performing the numerical calculation, we would 
like to point out that all the relevant energies such as the 
AF exchange energy, the kinetic energy and the energy 
for the formation of spin gap are contained in Eq.(^. 
Comparing the AF exchange energy in Eq.(^ to that in 
Eq.(|l]), we find two differences, a) there is an additional 
integral over v in Eq.(|^) and b) the temperature factors 
coth(ijj/2r) and 1/[1 — exp(— tj/T)] are different. At zero 
temperature or if the main contribution to the integral 
over to is dominated by those satisfying lu/T ^ 1, the 
temperature factors are identical. When v ~ \, Eq.(|7[) is 
the same as the renormalized spin susceptibility Eq.(^. 
Therefore, if the integral over v is dominated by that part 
around v k 1 and the integral over w is dominated by the 
spectral weight around w/T ^ 1, Eq.(^) can be reduced 
to Eq.(|l|). We will discuss this issue in detail in the first 
part of the following section. 
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III. NUMERICAL CALCULATIONS 

In this section, we carry out the calculations of the 
condensation energy based on Eq.(H) and Eq.(0). The 
purpose of this section is twofold. First, we will show 
that when a resonance peak occurs in the spin excitation 
spectrum, the result obtained by Scalapino and White 
Eq. (|l]) is almost identical to the condensation energy cal- 
culated from the perturbation diagrams shown in Fig. 2. 
Second, we will investigate the temperature and doping 
dependences of the condensation energy based on Eq. , 
and compare them with the experiments. 

Numerical calculations are performed by dividing the 
Brillouin zone into 128 x 128 lattices with t — 2J and t' = 
— 0.45t as usual. A quasi-particle damping F — 0.02J is 
used. We find that the results for the condensation en- 
ergy (including its magnitudes) nearly does not change 
when F is changed around 0.02J, though the spin sus- 
ceptibility indeed follows the change of F. The reason is 
because the change of F induces nearly identical changes 
for the spin susceptibility in both the normal and the SC 
states. 

The temperature dependences of the difference in the 
AF exchange energy calculated with Eq. (0) are presented 
in Fig. 3 for doping densities (5=0.16, 0.18 and 0.20, re- 
spectively. The results for the first term of the right 
hand side of Eq. (|^) with the same parameters are shown 
in Fig. 4. In all cases, both results are only weakly T- 
dependent at low temperatures, but fall off rapidly as 
T is increased towards T = 0.08J. Comparing Figs. 3 
and 4, one can see that the overall temperature depen- 
dences for both cases are remarkably similar for the dop- 
ing range investigated. Meanwhile, both have a simi- 
lar magnitude, though the latter is slightly smaller than 
the former. Thus, we may conclude that the first term 
of Eq.(||) and Eq.(^ may describe the same physics. 
Therefore, it provides a thermodynamical extension to 
the work by Scalapino and White As noted above, 
the most possible reason for the same temperature de- 
pendence is that the integral over v is dominated by that 
part around 1^=1 and the integral over frequency is 
dominated by that part satisfying cj/T 1. To see the 
variation of the integrand in Eq.(^) with v, we show in 
Fig. 5 the results for Imx"' {v,Q,uj), which is given by 
Eq.(0), at temperature T = O.OOSt and doping (5 = 0.16 
for several values of ly (1, 0.9 and 0.5). When ~ 1, the 
integrand x"' {v,q,uj) is identical to the spin suscepti- 
bility Eq.(^. The most remarkable change in the spin 
susceptibility across the SC transition is the emergence 
of the resonance peak at momentum Q = (tt, tt) below Tc, 
which is shown as the solid line in Fig. 5. In the frame- 
work of the d-wave BCS theory, the origin of the neu- 
tron resonance peak has been attributed to a collective 

p^ , p5t , which corresponds to 
{Q,uj) — and Ibixq^ (Q, t^) approaches 



spin excitation mode |12 



to zero. It is caused by the step-like rise of the imagi- 
nary part of the unperturbed spin susceptibility Imxo at 
its threshold as shown in the inset of Fig. 5, where the 
solid line denotes its real part, the dashed line its imagi- 
nary part and the dotted line represents —l/uriJjQ with 
1^=1. According to the Kramers-Kroenig relation, a 
logarithmic singularity in its real part Rexo occurs due 
to this step-like rise, which exhibits as a peak in the nu- 
merical calculations on finite lattices. This enhancement 
in Rexo shifts downward the position of the resonance 
mode, which is the cross between the dotted line and the 
solid line in the inset of Fig. 5, into the SC gap where no 
damping for spin excitations is expected. When v de- 
creases, the dotted line in the inset of Fig. 5 will rise, so 
the position of the cross will move to higher frequencies 
above the threshold where the dampings become nonzero. 
As v decreases further, there is no cross anymore. How- 
ever, if 1 + rjJjQUexoiQ,!^) is small, a peak still oc- 
curs. But the intensity of the peak will decrease, such 
as the dashed line for = 0.9 in Fig. 5. This may be 
regarded as the proximity effect of the resonance mode. 
As ly moves far away from 1, no peak exists and only a 
broad hump appears such as the dotted line in Fig. 5 for 
ly — 0.5. From the figure, one can see that the integrated 
spectral weight over uj in the SC state decreases with the 
decrease of v. However, little change occurs in its normal 
state results (the extrapolated value). Therefore, it is the 
existence of the resonance peak that the integral over v 
is dominated by that part around ly — I. Because the 
resonance peak is around uj — 0.52J( 0.53J and 0.52J 
for i5 = 0.18 and S = 0.20, respectively), so the tem- 
perature factors coth(a;/r) and 1/[1 — exp(— w/T)] are 
essentially equal to 1 even for the highest temperature 
T = 0.08 J considered here. As for the momentum distri- 
bution, we wish to point out two facts: a) the resonance 
peak appears at Q = (tt, tt) and drops rapidly as q moves 
away from Q; b) the summation over q is weighted by a 
factor 7g — cos{qx) + cos{qy) being maximum at Q. Ob- 
viously, both help to filter out the contributions around 
the resonance peak. Therefore, the very similar temper- 
ature dependence of the results calculated with Eq.(|^) 
and with the first term of Eq.(||) suggests that a large 
part of the condensation energy due to the AF exchange 
interactions comes from the spectral weight around the 
resonance mode, as first argued by Demler and Zhang 

Because the integral over ly in Eq.(H) is dominated by 
the contributions around v — 1 and only in this range of ly 
is the integrand Im^"' [v, Q,^^) approximately identical 
to the spin susceptibility Eq.(^), it leads the magnitudes 
of the condensation energy calculated with the first term 
of Eq. (|^) to be smaller than those calculated with Eq. (^ , 
as one can see from a comparison of Figs. 3 and 4. In 
Fig. 4, we also show the temperature dependences of the 
resonance peak intensities for doping S =0.16(squares), 
0.18(circles) and 0.20(triangles). For all doping levels 
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from the optimal to overdoped regimes, we find that the 
peak intensity and the condensation energy due to the 
AF exchange interactions follow the same temperature 
dependence. This result is in good agreement with the 
recent experiment on YBa2Cu3 06.93 which is an op- 
timally doped system, though the result for the over- 
doped system waits for future experiment testing. This 
coincidence further supports the above conclusion that 
the contribution to the condensation energy from the AF 
exchange corrections arises mainly from spin excitations 
around the resonance peak. 

The total condensation energy based on the renormal- 
ized RPA approach is plotted as a function of temper- 
ature in Fig. 6 for different dopings S = 0.16,0.18 and 
0.20. From a comparison with Fig. 3, we find that the 
results calculated here exhibit different temperature de- 
pendence, i.e., they decrease more rapidly with the rise 
of temperatures. This difference is due to the rapid 
variation with temperature of the contribution from the 
mean-field part. Meanwhile, their magnitudes are only 
about half of those calculated with Eq.(|^). We note that 
the condensation energy for doping S = 0.16(nearly op- 
timally doped) at zero temperature is around 0.023J, 
which is in reasonably agreement with the estimated 
value 0.016J based on the experimental data for the 
spin structure factor and the g-width of the resonance 
peak in YBa2Cu307 ||^ and also with the experimental 
data Ec = 6K(0.23 J « 30K if J = 1300K is taken) in the 
specific-heat measurement pO[ |. Thus, our investigation 
based on the renormalized RPA approach to the t — t' — J 
model gives a quantitative account for the condensation 
energy. 

The doping dependence of the condensation energy cal- 
culated with Eq.(H) is shown in Fig. 7. An obvious feature 
is that the condensation energy decreases as the dop- 
ing density rises from the optimally doped to the over- 
doped regime. This is related to the doping dependence 
of the maximum SC gap Aq at the mean-field level which 
shows the same trend, for example Ao=0.335J, 0.322J 
and 0.311J for 5 = 0.16,0.18 and 0.2 at T = 0.005 J, re- 
spectively. Moreover, an approximate linear variation of 
the condensation energy with doping is found. This re- 
sult is consistent with the experimental data determined 
from the heat capacity measurements |2^ . However, we 
note that the condensation energies shown in Fig. 6 do 
not approach to zero when the temperature is as high 
as 0.08J which may correspond to 104K if we choose 
J = 1300K. This is the fault of the slave-boson mean- 
field calculation which over-estimates Tc so that Aq does 
not approach to zero even at T = 0.08 J. Therefore, even 
though our approach gives good quantitative results near 
zero temperature, it may break down at temperatures 
near true transition temperatures. 



IV. CONCLUDING REMARKS 

In summary, based on the slave-boson theory of the 
two-dimensional t~t' — J model, we calculate the super- 
conducting condensation energy for optimally doped and 
overdoped high-Tc cuprates at finite temperatures using 
a renormalized random phase approximation. It is com- 
posed of the contributions from the mean-field Hamilto- 
nian and that from the antiferromagnetic spin fluctua- 
tion. 

We show that the contribution from the antiferromag- 
netic spin fluctuation is essentially the same as the ar- 
gument by Scalapino and White which relates the 
condensation energy to the difference of the antiferro- 
magnetic exchange energy between the normal state and 
the superconducting state, due to the existance of the 
neutron resonance mode. However, both the tempera- 
ture dependence and the magnitudes of the total con- 
densation energy after including the mean-field part are 
different from those due to the AF exchange energy. Our 
result for the total condensation energy gives a quantita- 
tive account for the observed condensation energy near 
zero temperature. We also calculate the doping depen- 
dence of the condensation energy and find it in agreement 
with experiments. 

Taking into account of our previous studies of the re- 
lationship of the resonance neutron peak to the angle- 
resolved photoemission spectra jl^ and the tunneling 
spectra [0, we believe that the resonant spin collective 
mode plays an important role in determining many phys- 
ical properties in the superconducting state of high- Tc 
cuprates. A particular interesting conclusion from this 
work is that the renormalized random phase approxima- 
tion seems to have already captured the main features 
of the resonance mode. However, it is introduced phe- 
nomenologically and its derivation is left for future work. 



ACKNOWLEDGMENTS 

JXL was supported by the National Nature Science 
Foundation of China. CYM thanks the support from 
NSC of Taiwan under Grant No. NSC 89-21 12-M-007- 
091 and TKL for NSC 89-2112-M-001-103. CDC thanks 
the supp ort by the Minist ry of Science and Technology 
of China (|nkbrsf-gl999064| 6). 



[1] G.V. Chester, Phys. Rev. 103, 1693 (1956). 
[2] A. Leggett, J. Phys. Chem. Sohds 59, 1729 (1998); Phys. 
Rev. Lett. 83, 392 (1999). 



5 



J.E. Hirsch and F. Marsiglio, |:ond-mat/9908322| (unpub- 
lished) . 

D. J. Scalapino and S.R. White, Phys.Rev.B 58, 
8222(1998). 

E. Demler and S.C. Zhang, Nature 396, 733 (1998). 

A. Abanov and A.V. Chubukov, Phys. Rev.B 62, R787 
(2000). 

J. Rossat-Mignod et al, Physica C 185-189, 86 (1991). 
H.A. Mook et al, Phys. Rev. Lett. 70, 3490 (1993). 
H.F. Fong, et al., Phys. Rev. B 54, 6708 (1996); Phys. 
Rev. Lett. 75, 316 (1995). 

J.W. Loram, K.A. Mirza, and P.A. Freeman, Physica 
C171, 243 (1990); 

J.W. Loram, K.A. Mirza, J.R.Cooper, W.Y. Liang, and 

J.M. Wade, J. Supercond. 7, 243 (1994). 

P. Dai, H.A. Mook, S.M. Hayden, G. Aepph, T.G. Per- 

ring, R.D. Hunt and F. Dogan, Science 284, 1344 (1999). 

J. Brinckmann and P.A. Lee, Phys. Rev. Lett. 82, 2915 

(1999). 

J.X. Li, C.Y. Mou, and T.K. Lee, Phys. Rev. B 62, 
640(2000). 

This was also the regime suggested by Scalapino and 
White 0, that could be used to check the validity of 
Eq.(§). 

G. Rickayzen, in Superconductivity, edited by R.D. Parks, 
Vol.1 (Marcel Dekker, NY) (1969). 

C.Y. Mou, J.X. Li, and T.K. Lee, in preparation. 

C.L. Wu, C.Y. Mou, and D. Chang, jcond-mat/0012168 
(unpublished). 

For a recent review, see T. Timusk and B. Statt, Rep. 
Prog. Phys. 62, 61 (1999). 

M.U.Ubbens and P.A.Lee, Phys. Rev. B 46, 8434 (1992). 
T. Tanamoto, K. Kuboki and H. Fukuyama, J. Phys. 
Soc. Jpn. 60, 3072 (1991); T.Tanamoto, H.Kohno and 

H. Fukuyama, tbid. 61, 1886 (1992). 

P.A. Lee, Physica C 317-318, 194 (1999). 

G.D. Mahan, Many-Particle Physics Chapter 3 (Plenum 

press, NY) (1990). 

J. R. Schriffer, Theory of Superconductivity, 
(W.A.Benjamin, Massachusetts, 1964). 
G. Blumberg, B.P. Stojkovic and M.V. Klein, Phys. Rev. 
B 52, 15 741 (1995); D.Z. liu, Y. Zha and K. Levin, Phys. 
Rev. Lett. 75, 4130 (1995). 

J.X. Li, W.G. Yin, and CD. Gong, Phys. Rev. B 58, 
2895 (1998). 

J.L. Tallon, J.W. Loram, G.V.M. Williams, J.R. Cooper, 

I. R. Fisher, J.D. Johnson, M.P. Staines, and C. Bern- 
hard, Phys. Stat. Sol.(b) 215, 531 (1999). 



FIGURE CAPTIONS 

Fig.l Fcynman diagram for the renormalized spin 
susceptibility calculated in the random phase approx- 
imation(RPA). The solid lines represent the fermionic 
Green's functions and the dashed line the antiferromag- 
netic coupling constant Jjq. 

Fig. 2 Feynman diagram for the contribution to the 
thermodynamic potential from antiferromagnetic fluctu- 
ations. The solid lines represent the fermionic Green's 
functions and the dashed line the antiferromagnetic cou- 
pling constant Jjq. 

Fig. 3 Temperature dependences of the difference in an- 
tiferromagnetic exchange energy between the normal and 
the superconducting states Eq. (|^) for hole concentrations 
6 =0.16, 0.18 and 0.20, which is argued by Scalapino and 
white Q to be the condensation energy. 

Fig. 4 Temperature dependences of the condensation 
energy due to the antiferromagnetic fluctuation correc- 
tion [the flrst term in the right hand of Eq.(9)] for hole 
concentrations S =0.16, 0.18 and 0.20. The solid symbols 
denote the intensities of the resonance peakes at different 
temperatures (squares for S = 0.16, circles for S = 0.18 
and triangles for 6 = 0.20). The intensity is scaled in 
order to compare it with the condensation energy. 

Fig. 5 Frequency dependences of Im^^' (i^, Q,Ci;) at 
temperature T = 0.005 J, Q = (tt, tt) and doping S — 0.16 
for 1/ =1, 0.9 and 0.5, respectively. Inset shows the bare 
spin susceptibility Imxo(Qi The solid line denotes its 
real part and the dashed line its imaginary part. The 
dotted line represents — l/77J7Q(see text). Similar re- 
sults are obtained for other dopings expected that the 
positions of the resonance peak are shifted, so are not 
shown here. 

Fig. 6 Temperature dependences of the total condensa- 
tion energy calculated with Eq. (^ for hole concentrations 
S =0.16, 0.18 and 0.20. They show more rapid drops in 
comparison to Fig. 3. 

Fig. 7 Doping dependence of the total condensation en- 
ergy calculated with Eq.(|) at T = 0.005 J. 
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Fig. 3, Li etal. 
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Fig. 6 Li, etal. 
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